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1. INTRODUCTION 


In this memorandum we discuss the algebraic connection between a 
rational transfer matrix H ( A ) and its realization by a time-invariant 
dynamic system, or triple of constant matrices (C,A,B). It is shown 
that both H ( A ) and its matrix realizations can be represented by a 
suitable morphism of modules over the ring of polynomials in X. Mor- 
phisms of this type are the fundamental algebraic objects in the theory 
of finite, constant linear systems. 

The discussion to follow was stimulated by results of Kalman, 
Gilbert and Youla (documented in Section 8) and makes no important 
claim to originality. However, our development differs in some res- 
pects from that of the authors cited and is presented as a basis for 
research on certain problems of multivariable system synthesis (e.g., 
decoupling) which have not so far been treated from a module-theoretic 
point of view. 

As algebraic prerequisite we draw only on the elementary theory 
of modules as given, for example, by MacLane and Birkhoff [11]. 

2. NOTATION 


R n is real n-dimensional vector space; R=r 1; R [ A ] is the ring of 
polynomials over R in a single transcendental element X; R n [X] is the 
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free RfXJ -module generated by R n . All modules which arise have Rf XI 
as associated ring. Dual modules are denoted by a prime, e.g. , (R n (X) ) ' . 
In general, script light-face will denote vector spaces over Rj script 
bold-face, modules over R[X] • Morph isms of R-vector spaces will be 
called maps and denoted by A,B, • • . Morphisms of modules are printed in 
bold-face S,T, ... A map BsRf^R* 1 has a natural extension to a morphism 
of free modules BiR^tXJ^R 11 !!] as follows! If {u, ,iem> is a basis for 
R* then the u^ also serve as a basis (set of independent free generators) * 
for R?*tX] , and we define Bu^Bu^, iem, regarding R n as a subset of 
R**(X]. Finally, if a,9eR[X], [0) a is the residue class of 9 (mod a) in 
che quotient ring R[X]/a(X)R(X] . 


3 . TRANSFER MATRICES AND MORPHI SMS 


A transfer matrix (TX) is any (finite) rational matrix of the form 
H(X) - a<X) -1 H (X) (3.1) 

A A 

where a(X)eR[X] and H(X) is a polynomial matrix with deg H< deg a. H(X) 
reduced if the GCD of a(X) and the elements of H(X) is 1. Every 
p*m TX determines a morphism of modules as follows* Let (Uj,jem) be a 
basis for R 01 and (y^/icg) a basis for R p . The elements 

y^X) * + o(X)R p (X), iep 

form an independent set of generators for the torsion module R P [X]/ 
a <X)R p [XJ. Define 


Hsrf^EXj R P CX)/a(X)R p [X) 


(3.2) 


by the rule 



P 

E : ., (X)y i (X), jcm 


(3.3) 


where H^fX) is the (i,j)th element of a(X)H(X). A morphism of R(X)- 
modules of type (3.2) will be called a transfer morphism (TM) . 
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Every TM determines a TX as follows. With generators , y^(A) 

chosen as before, the evaluations (3.3) of the Hu. determine polynomial: 
a * 3 

H. .(X) uniquely, mod a(A). In particular let H. .(A) be the unique mem- 

1J A lj 

ber of its residue class with deg < deg ot, and let H(A) be the matr: 
of H^j (A) . Then H(A) = ot(A) ^ H(A) is the TX determined by H. 

Let H be a TM as in (3.2). Clearly 

Ker H } a ( A ) R m [ A ] 

Let ]_ be the subset of polynomials BeR[A] such that 
Ker H D B (A)R m [A] 

Clearly I_ is an ideal of R [ A ] , so that 7_ = a (A)R[A] for a polynomial 

a Q eR[A] . In general a Q | ot . Taking a and a Q to be monic, we say that 

the TM H is reduced if a = a. It is easily seen that H is a reduced 
— o — 

TM if and only if the TX determined by H is reduced. Obviously, each 
TX is equal (as a rational matrix) to a unique reduced TX with monic 
denominator polynomial. To avoid ambiguity we regard a p*m TX as a 

A 

pair (a(A),H(A)) with a (A) monic; then the corresponding TM is uniquely 
determined after choice of bases in R* 3 and R m . 

4. DYNAMIC SYSTEMS 


For our purposes a dynamic system (DS) is a triple of maps (C,A,B) 


C ; R - R^, A : R - R 11 , B:R m - R 1 


(4.1) 


which can be thought of as coming from the equations x = Ax + Bu, y = 
Cx. Every DS determines a TX by the rule 


H ( A ) = C ( A-A) _1 B 


(4.2) 


where we have chosen bases for R m ,R n ,R^ and used the same symbols for 


- 3 - 


the corresponding matrices. 


Let T:R n -*■ R n be an automorphism of R n . Relative to fixed bases 
in R m ,R n ,R P the DS (CT,T -1 AT,T -1 B) determines the same TX as (C,A,B) . 

DS which are related in this way will be called equivalent (— ) . 

With C,A,B as in (4.1), let C,A,B be their natural extensions to 
morphisms of free modules 

C :R n [X j - R P [X] , A:R n [X] - R n [X], B:R m [X] -> R n [X] 

If a(X)eR[X] is divisible by the minimal polynomial of A, the morphism 

a (X ) (X-A) ” 1 : R n [ X ] - R n [X] 

is defined in the natural way. Let be the canonical projection 
: R p [ X ] - R P ( X ] / a ( X ) R P [ X ] 

We say that (C,A,Bj realizes the TM H if a(X) is divisible by the mini- 
mal polynomial of A, and 

H = P^Cct(X) (X-A) -1 B (4.3) 

It is clear that (C,A,B) realizes H if and only if the TX of H satisfies 
(4.2); and that equivalent DS realize the same TM. There is defined in 
his way a function from the quotient class {DS}/~ to the class of TM. 

5. CONSTRUCTION OF A CANONICAL REALIZATION 


Write 

U = R m [X] , y = R P [ X ] /a (X ) R P [ X ] 


and let H :U-*y be a TM. Let 
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(5.1) 


M = U/Ker H 


Lemma 4 . 1 

There is a unique integer n, and a map A:R n -► R n , such that 


R n [X]/ (X-A) R n [X] « M 


(5.2) 


Proof 

Let u(X) = u(X)+Ker H e M. Since cx (X ) R m [X] C Ker H, 

a(X)u(X) = a(X)u(X) + Ker H = Ker H = 0; 

so that M is a torsion module of finite type over the principal ideal 
domain R [ X ] . By known structure theorems ([11] ,Ch. 10, Ths . 5,6) there 
exist unique monic polynomials a^,...,a^ such that a^|a, a 2 | , . . . , 

a kl a k-l' deg a k ” 1 ' and 
k 

M » 0 R[X]/ot i (X)R[X] (5.3) 

i=l 


n = ^ deg 
i=l 

By the rational canonical structure theorem ([11], Ch. 10, Th. 8) there 
is a map A : R n -*■ R n , unique up to similarity equivalence, such that 

k 

R n [ X ] / ( X -A ) R n [ X ] s 0 R[X]/a.(X)R[X] ; (5.4) 

i=l 

the a^(X) are the invariant polynomials of A. By (5.3) and (5.4) the 
result follows. I 


We shall say tnat a realization (C,A,B) of H is canonical if 


(5.2) holds. 

It is convenient in the following to define M concretely as an 
isomorphic copy of the module on the right side of (5.4), namely 
k 

M - 0 B i (X)R[X)/a(X)R[X] (5.5) 

i=l 

where B = a/ou , iek. H can now be exhibited as a product 

H = Hq-I" (5.6) 


where Hj is the canonical projection of U on M and is the induced 
monomorphism from M to V , as in the diagrams below. 



(Dl) 




U ■ 


•M 


(D2) 


0 M ~V (D3 ) 

The subscripts 1,0 refer to "input" and "output". D2 and D3 are exact 
sequences, i.e., Hj is an epimorphism and is a monomorphism. It 
will be seen in Section 6 that "epic" and "monic" have the system — 
theoretic meaning of "controllable" and "observable". 

Write 


X - R n [ X ] / ( X -A ) R n [ X ] 


(5.7) 
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and A = (X-A)R n [X]. It is easy to check that elements xeX have the 
representation 

x = x + A, (5.8) 

where xtR n is uniquely determined by x. Furthermore, if x(X) = x + 
(X-A)x*(X) for xeR n , x*(X)eR n (X), then 

Xx (X ) = (A+ (X-A) ) x (X ) e Ax+ (X-A)R n [X] , 

so that 

Xx = Ax + A (5.9) 

Next, since Xa;M, '5.3) implies 
k 

X*0X i (5.10) 

i=l 

where X^ is a cyclic module with generator 

?i ■ 9i + A (5.11) 

and g^ has minimal annihilating polynomial ou(X). 

We can now exhibit explicitly an isomorphism T:M-*X. Write 
= 6 i (X)R(X]/ot(X)R[X] , iek 

and let e. be the generator of M. , i.e., e. = [8;] a . Define 

— i s — i *i ; e i ,— ie - (5.12) 

and finally 
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1 - ® Si 

i=l 

The situation is now summarized by the diagram: 


s x 


So 


u © ... © ^ -y 


Si 


Ik 


(D4) 


Si®. • .©^ 


The remainder of the realization procedure amounts simply to 
filling in the diagonals so that D4 commutes: 


Si 


So 


\ 

X 1 
J - \ 


\ 


s 

s' K 


(D5) 


Define J:U -*■ X as the morphism 


J = T Hj 


(5.13) 


Let {u^,iem} be a basis for R m , hence of £. Then 


JUi = T(H lUi ) 


for uniquely determined vectors b^eR n . Thus J determines a unique map 
B:R m *■ R n by the rule 

Bu i = b i' ^ e - (5.14) 
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The morphism K:X -*• V is defined by KT = or 

K » HT" 1 
— — o— 

= Mo t?.! 1 ©* • *© T k 1) (5.15) 

Observe that tT 1 maps each element x(X)eX^ into its corresponding re- 
sidue class in H.; i.e., if x = S(A)g^+A e X^, then 

T“ X x = [C (A) ] a (5.16) 

To determine an explicit expression for the T^ we make a digres- 
sion on cyclic maps. Let G:R n -*- R n be cyclic, with minimal polynomial 
Y ( A ) and generator geR n . Introduce the element h'e(R n )' conjugate to' 
g, defined by 

h'g = . . . = h'G n " 2 g = 0; h'G n_1 g = 1 <5.17) 


Let 

Y(X) = X n - (c.+c 0 X+ . . ,+c X^ 1 ) 
l t n 

for real numbers c^; and define polynomials 

Y (i >(X) = \ n ~ L - (c i+1 +c i+2 X+...+c n X n ~ 1 ' 1 ), ien (5.18) 

so that y (l) (X) = Xy ( 1+1) (X) - c i+1 , i£n-l, and y {n) (X) = 1. For the 
natural extension G:R n [X] ■* R n (X) it is easy to verify 

n 

Y ( X ) = (X-G) ^ Y U) (X)G i_1 
i=l 


and we define 
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(5.19) 


Y(X) (A-G)" 1 = ^ Y (l) (MG 1 ” 1 
i=l 


With h' the natural insertion of h' in (R n [A])', define the functional 
F:R n [X] - R [ X ] by 

F = h ' y (X ) (X-G)" 1 (5.20) 


Lemma 5 . 2 


For every £eR[A] 

F4 (G)g = C (X) mod y (X) 


(5.21) 


ys a 

If Fe (R n [ X ] ) ' has the property (5.21) then F-F = y (X)F for some 
F e (R n [ X ] ) ' . 

Proof 

To verify (5.21) note from (5.17), (5.19) and (5.20), 

Ty i] (X ) h *G 1 ~ 1 g = Y (r ‘U) = 1? 

i=l 

F[XG 1 -(X-G)G 1 ]g 
XFG 1 g - a(X)h'G 1 g 
XFG 1 g mod y (X ) ; 

and (5.21) follows by induction on i. The uniqueness (mod y) of the 
functional determined by (5.21) is implied by the fact that {G 1_1 g,ien} 
is a basis for R n [A].l 

Now let R n be the subspace 


Fg = 


„.i+l 
FA g = 
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X^ = {x:xeR n ,x+A e X^} 


Then X^ is A-cyclic with minimal polynomial a^(X) and (by (5 
tor g^ . Let X^[X) be the free module generated by X^. Then 
phism 

a.(X) (X-A)" 1 : X ± ( X ] - R n [X] 


is well-defined, 
by 


y«i ■ 


Next, write n^=dim X^ (=deg cu) and define 






1 


Finally, let F. :X^[X1 ■+ R[X] be the functional given by 
F\x(X) = h|a i (X) (X-A) -1 x (X ) , x(X)eX i [X] 


Lemma 5 . 3 


,-l 


The morphism T^ : X^-*-M^ is given by 


T? 1 (x) = [6 i (X)F i x] a 


for x = x + A e X . . 

— —l 


Proof 


Since xeXX X^[X], (5.23) makes sense. Now apply Lemma 
X i in place of R n , the restriction of A to in place of G, 
a i (X) . | 

It remains to compute K from (5.15) and (5.23). Let {y 
a basis for R p , so that 


y i (X) = y.^ + a(X)R p [X], iep 


11)) genera 
the mor- 


•i e (— i [X] ) ' 


(5.22) 


(5.23) 


5.2, with 
and y (X ) = 

. ,iep} be 
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is an independent set of generators for the torsion module V_. Applying 

H to the generators e. = [B.,] of M, there results 
—o j D ® 

p 

^0*3 = j£ - (5-24) 

i=l 


for unique (mod a) polynomials n\ j . It will be shown that In 

fact with aj=a/Bj, we have ot^e^O e M^, so that 


0 = S 0 (Oj e jJ 


- “jSo'V 


p 

" a j (X) Z/ r ij U)y i 

i=l 


(X) 


Since the y^(X) are independent generators, (X) ri\j (X) = 0 (mod a(X)), 
i.e. , 


a j(X)^ij ( X) e a (X ) R l X ] 

°r j ( X ) e Bj(X)R[X], and the assertion follows; so 

P 

Mo e j = B j (X) 2 n ij (X)y i (X) ' je - (5 * 25) 

i=l 


for suitable n^eRfX]. Also t by (5.21), 

h j n i j Oj (X) ( X —A) X g j = n i j(X), mod , ie£, jek (5.26) 

k 

Now let x+AeX, so that x = ^ x. with x . eX . , jek, uniquely deter- 

j=l - 1 3 ^ 

mined by x. Let {y|,ie£} be the basis dual to {y.,ie£> in (R p ) 1 ; define 
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(5.27) 


the map C:R n -*• R P by the rule 

k 

y i Cx = ^^ h j n ij (A ) x j • ie E' xeRn ; 

j-1 

and let y|e(R P [X])', C:R n [X] ■+ R P [X] be the corresponding natural ex- 
tensions. Then 


K (x+A) = H (T -1 (x+A) ) 
— — — o — — 


k 

= H T^[h!ct. (X) U-A) _1 x.]e 
-o Z-r -3 3 ~ 3 ! 

i=l 


k p 

= ^[h! aj (X) (X-A)“ 1 x j ie j ( x )^ n ij 
j=l 


i=l 


(by (5.25)) 


P k 


EZ [htriij (A)B. (X)a. (X) (X-A) _1 x . ] y £ (X ) 

i=l j=l 


(by (5.26).) 


t- — 

= ^[y! ^Ca(X) (X-A)" 1 x j ]y i (X) 


i=l j=l 


(by (5.27)) 


P 

= y^ly^Ca(X) (X-A) ~ 1 x]y i (X) 
i=l 


= P Ca (X) (X-A) _1 x 


(5.28) 


for all x+A e X. 


To 


complete the discussion let u^eR m be a free generator of U 


Then 
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(by (5.28)) 


Hu^ = K Ju. 

= K(Bu i +A) 

= P^Ca ( X ) (X-A)" 1 B Ui 
= P^Ca (X) (X-A) _1 Bu i , iem; 

so that 

H = P^Ca ( X ) (X-A) _1 B ; 

i.e., (C,A,B) is a realization of H. 

We have shown that a canonical realization of a TM always exists. 
It will now be shown to have important special properties. 

6. PROPERTIES OF CANONICAL REALIZATIONS 

Let (C ,A,B) be a DS. The pair (A,B) is controllable if 
n 

J^Range (A 1 “ 1 B) = R n ; (6,1) 

i=l 

and (C,A) is observable if 
n 

P|Ker (CA 1_1 ) = 0 * (6.2) 

i=l 

These well known concepts, due to Kalman, are fundamental. The DS 
(C , A , B) will be called complete if (A,B) is controllable and (C,A) is 
observable. It is easily checked that completeness is a class pro- 
perty relative to equivalence (~) . 

From a practical viewpoint it is of interest to realize a TM in 
such a way that the state space R n is of minimal dimension. A 
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realization with this property is nunimal . 

The central result of realization theory is 
Thoorom_ b_. 1 

Let the PS be a real ization of the TM H. The t ol lowing 

properties of (C , A, N) ajre equ ivalent: (i) ( C , A , B ) incanon leal ; (ii) 

(C,A,t’>) is complete; (iii) (C,A,B) is minimal, furthermore , therc^ is 
a bi jection between t ho class of reduced TM and the quotient class of 
complete PS modulo equivalence (*») . 

In the following, the notation of Section 5 will be used freely. 

Lemma_6_. 1 

Let S:\*\ be a morphism and define the function S:R n> R n by the 
> \ 

rule: S(x+A) = x + A implies Sx = x. Then S is a map and SA -- AS.' i_f 
S is an autom orp hism, so is S. 


rroof 


Since x, x t R n are uniquely determined by the elements x + A, x + 
A ;■ the function S is well-defined and is clearly a map. Recalling 
( 5 . 9 ) , 

S (Ax + A) = S ( \ (x + A) ) 

= \S (x + A) 

= ,\ (x + A) 

= Ax + A 

so that S(Ax) = Ax = A(Sx). The second statement of the lemma is 
obvious . I 

Lemma 6 . 2 

If (C,A,B) is a canonical realization, it is complete. 
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Proof 


Suppose first that (C,A,B) is constructed as in Section 5, i.e., 
B is given by (5.14) and C by (5.27). As for controllability, since 
Hj is epic there are polynomials (X) (iek,jem) such that 



m 


£ 

■4 = 1 



(X)u. 


m 

y; pij ( x )Hx ( u j ) ' iE ii 

j=i 


Therefore 


g . + A = Te . 

- — i 

m 

“ £ P ij (A)T(H I u j ) 

3=1 


- £ lOij CA>bj + AJ; 
J-i 


m 



j=l 


Since any xeR n can be written 
k 

x = S C i (A)g i 

i=l 

for suitable £^eR[X], controllability is proved. 

For observability, since is monic and T is an isomorphism, 

- 1 ~ ° 

HqT is monic, so that 

x = x + A e X and H^T ^ x = 0 imply x = 0 (6.3) 

By (5.28) and (6.3), ^ x=0 if and only if the polynomials 
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mod a (X ) , ie£ 


y 'Ca (A ) (A -A) ^x = 0 , 
i 

or 

Ca (A ) (A-A)" 1 x = 0 (6.4) 

in R^[A] . By use of (5.18) and (5.19) it is seen that (6.4) is true 
if and only if. CA 1 1 'x=0, icn; that is, 

xeR n and CA* *x=0, ien, imply x=0, 
so that (C , A) is observable. 

/N A 

Now let (C , A , B) be any canonical realization, i.e., (5.2) holds 

for A. Define morphisms 

A A 

J:U+X; u h- Bu + A, U£R m (6.5) 

K:X-*y ; x + A K P^Ca(A) (A-A) -1 x, x + A e X (6.6) 

Then H=P^Ca ( A ) (A-A) ^B , so that 

Hu. = P Ca(A) ( A -A) -1 Bu . 

— 1 — Ot — — — 1 

= P^Ca ( A ) (A-A) _1 Bu i 
= K(Bu i +A) 

A A 

= KJu^ iem 


and therefore H=KJ . 

Next determine B,C according to (5.14) and (5.27), so that 
(C, A, B) is complete, and H=KJ, where J,K are given by (5.13) 
and (5.15). It will be shown that there exists an automorphism 
S:X+X such that diagram D6 commutes. Lot Kj : X+X/Ker K be the 
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a 

X 

I 

| - . (D6) 

X 

A A 

canonical projection, so that K = K 2 *Si' where K^:^_/Ker K -*> X is monic. 
Since K is also monic (by the construction of Section 5) , there exist 

A A A 

.submodules X n SX/Ker K and X z X such that X ZK z HU. Since J : U-*X is 
epic r a 1 HU - K ( JU) it follows that HU = KX, so that X^z X . As in the 
proof of Lemma 6.1, an isomorphism QrX^+X induces a map Q:X q ->X of the 
concur -Ting vector spaces; it is immediately verified that Q is an 
isomorphism, so that X q = X, hence = X, and X Z X/Ker K. If C is the 
vector subspace 

C = {x:xeR n , x+A^ e Ker K} 

/\ 

then, as before, an isomorphism Q:X+X/Ker K induces an isomorphism of 
vector spaces Q:R n -*R n /C, which implies C=0 , and so Ker K = 0. Thus K 
is monic, so there is an automorphism S of X such that KS '= K. Then 
KJ = h = K(S“lj); K monic implies J=S"^J; and D6 commutes, as asserted. 




By Lemma 5.1, S induces an automorphism S:R n +R n such that SA=AS 
Clearly B=S *B; and KS = K implies 


P^CSa (X ) (X-A) 1 = P Ca(X) (X-A) _1 S 


- P„Ca(X) (X-A)' 1 . 

— Ot — — 


By (5.28), 

y^CSa(X) (X-A) _1 x = ^Ca(X) (X-A) -1 X e R(X] 
for ie£ and xcR n . Hence equality holds for x ( A ) e R n [ A ] ; setting 
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x ( X ) ** (X-A)Xj(jen) for a basis {Xjijcn} of R n , there follows y^CS = 
y|C, ie£, and so CS = C. Thus (C,A,B) = (CS,S -1 AS,S _1 B) (C,A,B) and 

therefore (C,A,B) is complete. I 

Lemma 6 . 3 


Any two canonical realizations of a TM are equivalent. 


Proof 


If (C^/A^B^ and (C 2 ,A 2 ,B 2 ) are canonical realizations of H:U-*y, 
there are integers n^,n 2 such that 


n . n . 

R [ X ] / (X -A i ) R 1 [ X ] a U/Ker H 


By the rational canonical structure theorem n^=n 2 =n, say, and there is 
an R-vector space automorphism S of R n such that SA^A^S. Thus • 

(C 2 ,A 2 ,B 2 ) = ^c 2 ,s" 1 a 1 s,b 2 ^ ~ ^c 2 s,a 1 ,s" 1 b 2 ^ 

By the method of proof of Lemma 6.2, there are maps B,C such that 
^C 2 S,A 1 ,S“ 1 B 2 ^ - (C,A,B) 


and 


(Ci , A^ , B^ ) (C , A^ , B) ; 


and the assertion follows. I 


Lemma 6 . 4 

If (C , A , B) is a complete realization of a TM H, then it is canoni - 
cal . 

Proof 

In the notation of Section 5, it must be shown that 
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U/Ker H # X 

For this it is enough to exhibit an epimorphism S:U-*-X such that Ker S = 
For H. Define S by the rule Su i = Bu i + A, icm. Since (A,B) is con- 
trollable, S is epic. Let u(X)EKer S, so that Bu(X) = (X-A)u*(X) for 
some u*(X)ER n [X], Since (C,A,B) is a realization, 

Hu (X ) = P^Ca ( X ) ( X - A ) ~ 1 ( ( X - A ) u * ( X ) ) 

= (X ) a* (X ) 

= 0 

and therefore Ker SC Ker H. 

For the reverse inclusion, suppose 
P^CcMX) ( X-A) -1 x (X ) = 0 

for some x(X)eR n [X]. Since x (X ) =x+ (X-A) x* (X ) for some xeR n , there 
follows 

P^Ca ( X ) (X-A) -1 x = 0 
and therefore 

Cot (X ) (X-A)" 1 x = 0 

From this it follows by observability that x=0 , i.e., x (X ) e (X-A) R n [X] . 
Thus if u (X ) e Ker H then Bu(X) = 0, so that Su(X) = Bu(X) + A = 0 in X. 
Hence Ker HCKer S, which completes the proof. I 

Lemma 6 . 5 


A realization is minimal only if it is complete . 


-:o- 


Proof 


Let (C,A,B) be a realization of H, with A a map of R n . 

Write 

n 

E = S Range • 

i=l 

define A^jE+E, A^ = A | E ? and B^:R n -*E, uh*Bu. Choose f-CR n arbitrarily 
such that E0 F « R n ; let P:R n -*-E be the projection on E along F? and 
define C 1 :E-*-R p by the condition C^P = CP. We claim that (C^A^B^) 
realizes H. If u(A)eR m [A] then 

a (X) (A-A^) _1 BjUfA) = a (X ) (X-A) ” 1 Bu(A) e E[X] 

where £[A] is the free module generated by E. Since PE [A] C E.[A] 

Ca(X) (X-A) Bu (A ) = C^ (X ) (A-A^ " 1 BjUfA) 


and the assertion follows. 


From this we c nclude that (C,A,B) is a minimal realization only 
if E=R n , i.e., (A,B) is controllable. Similarly, by introducing as 

new state space the factor space 



(CA 


i-1 


) 


it can be shown that (C,A,B) is minimal only if (C,A) is observable. I 
Lemma 6.6 


A canonical realization is minimal. 


proof 


Any canonical realization has state space R n determined uniquely 
by the module 11 / Ker H. Now a minimal realization always exists (be- 
cause a realization does); by Lemma 6.5 it is complete; hence by Lemma 
6.4, canonical; and the result follows. I 

Proof of Theorem 6.1 


The first assertion follows immediately by Lemmas 6. 2-6. 6. For 
the second, if two complete DS realize the same TM then (Lemma 6.4) 
they are canonical, hence (Lemma 6.3) equivalent. On the other hand 
every TM has a canonical, hence complete realization. Finally, let 
(C ,A,B) be complete and realize the reduced TM H, i.e., 

H = P^Ca (A) (A -A) _1 B 


It will be shown that a (A) is uniquely determined as the minimal poly- 
nomial of A. Since (C,A,B) is canonical, H = KJ with K, J defined as 
in Section 5. Let be the minimal polynomial of A. If 

u (A ) = a x ( A ) v (A ) e a 1 (A)R m [A] 
m 

with v ( A ) = 6 i (A)u i , then 
i=l 


Ju(A) = (A) Jv (A ) 


= a x (A) 


£ e i (A)b i + a 
L i=l 


m 


= a x (A) ^ 9 i ( A ) b i + A 


i=l 


= 0 ; 
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therefore (A) R m [A] C Ker H. Since H is reduced, ala^; and by defini- 
tion of a realization, | ct . I 

7. COMPUTATION OF A CANONICAL REALIZATION 

We indicate briefly how the construction of Section 5 can be per- 
formed explicitly. Start with the p*m transfer matrix 

H (A ) = a ( A ) _1 H ( A ) 

A 

As a polynomial matrix, H(A) can be factored by elementary row and 
column operations, into the form 

H (A) = P (A) D ( A ) Q (A ) 


where P(p*p) and Q(m><m) are invertible polynomial matrices (i.e., have 
R-valued, nonzero determinant) and D(A) = diag [6^ (A) , . . . ,6 r (A) ,0, . . .0] 
is a p*m polynomial matrix with the diagonal shown (of length min(p,m)), 
and zeros elsewhere. The 6 ^ can be chosen such that 6 ^ | 63 » . . . 6 r _^ I • 
and are then unique up to factors in R. Then 

r 

H ( A ) = a ( A ) -1 ^ P i (A)6 i (A)q i (A) 
i=l 


■L 


P i (A)6 i (A)q i (A) 

aTTM 


^ Pi (A)q.(A) 

“ 2-r a. (A) + 

i=l 1 


E(A) 


^ V-J A A 

where p^,q^ are respectively the i column, row of P,Q; = a/GCD( 6 ^,a) 

6 . = 6 ./GCD (6 . , a) ; p.,q. are the remainders after division of 

^i^i' ^i' respectively, by et^; E is a polynomial matrix; and terms in 
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the summation have been relabelled to retain only 

(i.e. , k < r) . 


Suppose deg = n^. Let A^ be a cyclic n^'n^ 

form 

f 0 1 0 .... 0 j 
0 0 1 0 ... 0 I 


0 

a il 


0 1 

.a . 

in. 


where the a^ are given by 


n. 


»i«'> - ^ 


j=l 


(j) 


From a i (\) compute the polynomials a^ J (\)(jfn^) as 


fine n^-vectors 


( n i“A 

p| ( A ) = col 1 1 , \ , . . . , \ I 
qj(\) = row (X),...,^" 1 * (X)j 


Then it is easilv checked that 


-1 


a i (.\)(\-A i ) A = pt(\)qj(,\) + a i (X)E.(X) 


where E^(\) is a polynomial matrix. 

Next compute real matrices (n^m) and O (p*n 

CjPjU) = p ± ( \ ) , q* ( X ) B^ = q ± ( X ) (iek) ; 
Ci and B^ always exist. Finally, let 


of positive degree 


matrix in standard 


in (5.18), and de- 


(7.1) 


. ) such that 
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A = diag [A^ # • • • , Aj^] 


B, 


B = 


1 C — [C^ 1 • 1 • |Cj^] 


B, 


It is a straightforward exercise in matrix algebra to check that the 
foregoing construction mimics the abstract procedure of Section 5. In 
particular it is helpful to note that we can take g^ = col (0 , . . . ,0 , 1) 
and h| = row (1,0,... ,0) in each n.xn^ block. 


8. BIBLIOGRAPHICAL NOTE 


In [1] Kalman described the canonical structure of a dynamic sys- 
tem (relative to controllability and observability) which we exploit 
in proving Lemma 6.5. Gilbert gave in [2] a realization procedure for 
a transfer matrix having simple poles; the computation, via residue 
matrices, is straightforward. Then Kalman [3] observed that, in general, 
a minimal realization could be computed by starting with an arbitrary j 

realization and reducing to a complete realization, as in [1] ; the 
equivalence of minimality and completeness was announced but not proved. 

In [4] Kalman gave the first general prescription for construction of 
a minimal realization, identifying the state space dimension as the 
(McMillan) degree of the transfer matrix; a realization was computed 
from a prime-factor (i.e., Jordan) decomposition of the state space. 

The first published proof of the equivalence of minimality and complete- 
ness, and the fact that minimal realizations are equivalent (~) , is due 
to Youla [5] , who worked directly with the differential equation for 
the impulse matrix; the technique is more analytical than algebraic. 

The module-theoretic viewpoint appears in Kalman's fundamental paper 
16]; the approach is via an input-output description in the spirit of 
automata theory, with the states identified as Nerode equivalence 
classes. As shown in [6], there immediately result the modular structure 


of the state space# and the crucial factoring of the input-output map, 
of which our diagram D1 is a modified version. Suggestive remarks on 
the system-theoretic meaning of homology sequences can be found in the 
expository paper [7]. In [8] Kalman gave an alternative concrete reali- 
zation procedure by way of rational canonical decomposition of the state 
space? this method is summarized in our Section 7, where the formula 
(7.1) is the "main lemma" of [8]. We note that in [8] the term "canoni- 
cal" is used to mean what we call "complete". The algebraic relation 
between non-minimal realizations is discussed in [9]. Finally Ch. 10 
f [10] contains a systematic exposition of much of the work just cited, 
wOgether with a presentation of the B. L. Ho algorithm for minimal reali- 
zation of an impulse response matrix. 

In the present memorandum our approach follows [10] in taking as 
"input" module the free module R m [X] , but differs in the choice of 
"output" module as R P [ X ] /a (X ) R P [X ] rather than the free module 
of formal power series over R p (see definition (2.1) of [10] , Ch. 10) .' 
This set-up seems better matched to the "finiteness" of the transfer 
matrix as a purely algebraic object (with no particular interpretation 
attached to it as an input-output map) , builds torsion in right at the 
start, and perhaps results in a more symmetric role for the C matrix. 

The formal definition of a transfer morphism is new. As a final detail, 
the definition (5.17) of a "conjugate" generator, and the useful formula 
(5.21), do not seem to have appeared before explicitly. 
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